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Abstract 

This work is a theoretical investigation of the stability of the non-linear behavior of an oscillating 
tip-cantilever system used in dynamic force microscopy. Stability criterions are derived that may 
help to a better understanding of the instabilities that may appear in the dynamic modes. Tapping 
and NC-AFM, when the tip is close to a surface. A variational principle allows to get the temporal 
dependance of the equations of motion of the oscillator as a function of the non- linear coupling term. 
These equations are the basis for the analysis of the stability. One find that the branch associated to 
frequencies larger than the resonance is always stable whereas the branch associated to frequencies 
smaller than the resonance exhibits two stable domains and one unstable. This feature allows to 
re-interpret the instabilities appearing in Tapping mode and may help to understand the reason 
why the NC-AFM mode is stable. 
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I. INTRODUCTION 



During the last decade, Dynamic Force Microscopy (DFM) used in the Tapping mode has 
been found as a suitable tool to investigate surface morphology, particularly of soft materials 
and has been widely used to investigate a very large range of samples including polymersi*^, 
biological materials^ or organic layera^*^. In recent years, the DFM was developed around 
the Non-Contact Atomic Force Microscopy (NC-AFM) mode and has shown that contrasts 
at the atomic scale could be achieved on semiconductors and insulators surfaces^^^iiSiiiiiSii^. 
Experimental and theoretical features dedicated to the description of these two dynamic 
modes have been widely discussed in previous paper a^'^i^^i^^i^'^'i^^i-^^i^" . It was shown that the 
large sensitivity of the oscillating tip-cantilever system (OTCS) was based on the large value 
of its quality factor and on its non-linear dynamics in the vicinity of the surface^i*^^. Even if 
the precise origin of the NC-AFM images contrast remains an open question. Tapping and 
NC-AFM results are the consequence of the same non-linear behavior of the OTCS so that 
their differences are purely technical. The aim of this article is to show from a theoretical 
point of view that the non-linear dynamics of the OTCS leads to various stability domains 
that may help to a better understanding of the way the instabilities appear in Tapping and 
the reason why the NC-AFM mode, while being so sensitive, keeps, in most of cases, a stable 
behavior. 

The experimental study of the stability of the non-linear behavior of the OTCS can 
be performed through the study of approach-retract curves, from which are deduced the 
variations of the oscillations properties as a function of the tip-sample distanc o^'^i^^ . But 
the ultimate goal of the approach-retract curves is to extract the surface properties from 
amplitude and phase variations (Tapping) or from frequency shift and damping signal (NC- 
AFM). In other words, the measure of the evolution of the amplitude with a fixed drive 
frequency and a drive amplitude or the one of the frequency shift are two complementary 
ways to investigate the non-linear behavior of the OTCS. Thus in this work, the analysis of 
the stability of the oscillations will be done from the study of the distortion of the resonance 
peak of the oscillator due to the coupling with the surface. 

The paper is organized as follow. The first part is dedicated to a description of the 
non-linear behavior of the OTCS at the proximity of the surface. To do so, a specific 
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theoretical frame giving the exphcit temporal dependance of equations of motion of the 
OTCS is developed. The description is then based on the equations of the stationary state 
of the OTCS that can be interpreted either for the Tapping or for the NC-AFM. In the 

second part, these equations are used to analyze the stability of the stationary state. The 
third part is a discussion of the results obtained and the way they can be interpreted for 
DFM experiments. 



3 



II. NON-LINEAR BEHAVIOR OF THE OTCS AT THE PROXIMITY OF A SUR- 
FACE 



A. General features 

The sketch of the OTCS and notations used in the present work are shown in figffl When 
a tip oscillates above a surface, several situations can occur. If the oscillation amplitude A is 
smaller than the distance D between the surface and the equilibrium position of the OTCS, 
the tip never touches the surface. On the contrary, for oscillation amplitudes larger than 
D, the tip exhibits intermittent contact situations. The time during which the tip touches 
the surface depends on several factors, among which, the local stiffness of the surface is 
a key parameter. For example, for very soft materials, the indentation depth can be as 
much as the oscillation amplitude^i, thus the tip spends half of the oscillation period into 
the substrate. For hard materials, the time during which the tip touches the surface is 
negligible compared to the oscillation period. Thus, when intermittent contact situations 
occur, various assumptions are required to describe the contact between the tip and the 
surfac o^^i^^i^^ . For the present purpose, such an analysis introduces useless complications, 
therefore we only focus on the non-contact situations. When uniquely non-contact situations 
take place, the tip never touches the surface during the whole oscillation period so that the 
interaction between the tip and the surface can be described with a simple expression. 

The attractive coupling force between the tip and the sample is assumed to derive from 
a sphere-plane interaction involving the disperse part of the Van der Waals potential^''. The 
differential equation describing the position of the tip's function of time, z (t), is 

given by : 

* 

m*Z {t) + (t) + KZ (t) = J^exc COS (ujt) - yVint [z it)] , (1) 

with : 

Vint \Z it)] = , , (2) 

^nt[ \ )\ Q[D-Z (t)] ^ ' 

cuo, <5, m* and kc = m*u!Q are respectively the resonance pulsation, quality factor, effective 
mass and cantilever' stiffness of the OTCS. J-'exc and u are the external drive force and drive 
pulsation. H, R and D are the Hamaker constant, the tip's apex radius and the distance 
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between the surface and the equihbrium position of the OTCS (see figH)). 



Different ways can be used to observe the non-hnear behavior of the oscillating nanosphere 
at the vicinity of a surface. The most informative is to record several resonance curves as 
a function of the distance /jiLiSiSSiS^. Another possibility is to keep the oscillation at a 
given frequency with a fixed excitation amplitude and to record variations of the oscillation 
amplitude and of the phase as a function of D. This latter measurement is readily done by 
performing approach curves with the Tapping mod o^^i^^i^^i^^ . As noted in the introduction, 
the difference between Tapping and NC-AFM is purely technical but, from a practical point 
of view, a description of the variation of the oscillating behavior showing the correspondence 
between the two modes is not straightforward. For example, the NC-AFM mode is probably 
the simplest way to describe a DFM experimenlj^^, but the non-linear behavior appears in a 
quite subtle way while it is immediately observed with the Tapping mode. 

The present section is divided into three parts. The first one details the specific theoretical 
frame for the obtention of the equations of motion in amplitude and phase of the OTCS. 
A coarse-grained method gives the equation describing the time evolution of the stationary 
state of the OTCS as a function of the coupling term between the tip and the surface. This 
allows to compute the stability of the stationary state. The next part is a description of the 
distortion of the resonance peak as a function of the distance. This provides the basis of the 
discussion about the oscillating behavior and the stability of the branches which is detailed 
in section IIIII This also provides an easier way to discuss experimental measurements as 
those obtained in Tapping and NC-AFM modes. As an example, we use the evolution of the 
resonance peak to discuss typical variations of the amplitude and phase in Tapping mode 
when the drive frequency is slightly below the resonance one^. These variations, regularly 
observed, are the most obvious experimental evidences of the non-linear dynamics of the 
OTCS. The third part deals with the analysis of the resonance frequency shift observed in 
NC-AFM and the way its stability can be interpreted. 

B. Theoretical frame 

We search a solution to the temporal evolution of the OTCS by using a variational 
solution based on the principle of least action. Even though this approach exploits the same 
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physical concepts than the one which had led to the coupled equations in amplitude and 
phase of the stationary state of the OTC S^'^'i^^i^^ , it appears to be more general since here, 
the temporal dependance is explicitly obtained. We start from the definition of the action 
of the OTCS coupled to an interaction potential : 

S= [' C{z,z,t)dt, (3) 

Jta 

where C is the Lagrangian of the systen>2^ : 



C{z,z,t) = T-V + W 



(4) 



-kcZ {tf - Z (t) Texc COS iijjt) + Vint [z if)] 



m Uq 



Q 



z{t)z{t) 



Due to the large quality factor, we assume that a typical temporal solution is on the form : 

z{t) = A it) cos \ujt + if it)] , (5) 

where A (t) and (t) are assumed to be slowly varying functions with time compared to the 
period T = 2ti/uj. The underlined variables of z_{t) : 

z{t) = A it) cos [ut + ip_ it)] - A{t)[uj + ^ it)] sin [ut + ^ {t)] , 

are calculated along the physical path, thus they are not varied into the calculationa^S. 
While among the trial functions solutions that can be chosen the present one is not the 
more general, the principle of least action ensures that functionals of this type are the best. 

The equations of motion in amplitude and phase of the OTCS are obtained by considering 
the following coarse-grained operation. Let's assume a long duration At = t}, — ta with 
At ^ T and calculate the action as a sum of small pieces of duration T : 



S 



l-{n+l)T ^ / 1 \ „ 

J2 / Ciz,z,t)dt = J2\7f / Ciz,z,t)dt]T = J2^eT 

n „ y JnT J 



(6) 



Ce is the mean Lagrangian during one period and appears as an effective Lagrangian for 
a large time scale compared to the period. Owing to the quasi-stationary behavior of the 
amplitude and the phase over the period, the effective Lagrangian is calculated by keeping 
them constant during the integration. The calculations give : 



£e {A,A,(f,(p 



m 
T 



A^ + A\u + 



m ujq 



4 2 
^474 cos — AA (^uj + Lp^ sin — 



Vint [z {t)] dt 
(7) 



Note that the effective Lagrangian is now a function of the new generahzed coordinates A, ip 
and their associated generahzed velocities A, ip. At this point, remembering that the period 
is smah regardless to At = — ta during which the total action is evaluated, the continuous 
expression of the action is : 



S 



(8) 



where the measure dr is such that T <^ dr <^ ^t. 

Applying the principle of least action ^S* = to the functional C^i we obtain the Euler 
Lagrange equations for the effective Lagrangian : 

d ( dCp\ dCf. , d f dCe\ dC 

and — 







(9) 



dr \dA J OA dr \ dcp J dcp 

The amplitude and phase equations of motion of the OTCS coupled to an interaction po- 
tential Vint [z it)] are : 



A 



+ 



1 



^ _ ^ ^exc COS (V?) _ UJ dI{A,Lp) 

Q kc Tike OA 



2A 1 W ^ 



J^excsm{ip)l uj 1 dI{A,ip) 



(10) 



A irkc A^ dip 



with 



I{A,^) 



Vint [Z (t)] dt 



(11) 



The system given by equsllOlis defined for any kind of interaction potential and no particular 
hypothesis are required to perform the calculations. This implies that, provided that I {A, ip) 
is analytical, the equations of motion can be obtained for any kind of non-linearity. For 
instance, in ref.— , analogous calculations were performed for the Buffing's oscillator. If 
we consider the sphere-plane interaction involving the disperse part of the Van der Waals 
potential (see equEI), it was shown^'^ that : 



I{A,^) 



ttHR 



3a; _ ^2 



(12) 



Thus 



dl (A, if) TcHR 



A 



OA 



3a; (L)2_^2) 



3/2 



and 



dip 



(13) 



7 



Using usual dimensionless notations^^, the equsHUJof the coupled equations of motion be- 
comes : 

r. , .^2 -il a C0S(V9) aKa 

/2d IV sin(v;) ' ^-"^^ 

\a Q J aQ 

d = D / Aq the reduced distance between the location of the surface and the equilibrium 

position of the OTCS normalized to the resonance amplitude Aq = QJ-'exc/^c, a = A/Aq is 

the reduced amphtude, u = u / uq is the reduced drive frequency normalized to the resonance 

frequency of the free OTCS and = HR/ (^c^o) is the dimensionless parameter that 

characterizes the strength of the interaction. The explicit dependance of the coupling term Ka 

with the oscillation amplitude through the power law Aq^ indicates the non-linear character 

of the dynamics. 



C. Resonance peak and amplitude variation recorded during a Tapping experi- 
ment 

The equations of motion of the stationary solutions a and (p are obtained by setting 
(2 = = and d = = in equ|T31and lead to two coupled equations of the sine and 
of the phase of the OTCS previously calculated^^ : 



cosme 



cos {(f) = Qa{l — u^) 



aQua 



3 (d2 - a 



2^3/2 



(15) 



sin ((/?) = —ua 

Solving equHni gives the relationship between the sweep frequency and the amplitude at a 
given distance d^^ : 



u± (a) 



4Q2 



It 



\ 



1-4Q'' 



3(d2 



,2^3/2 



(16) 



The signs plus and minus are deduced from the sign of cos (</?) and correspond to values of 
the phase ranging from to —90° («_, cos (ip) > 0) or from —90° to —180° («+, cos {(f) < 0), 
in agreement with the sign convention of the phase in equ|Hl From equ^l is calculated the 
resonance peak at any reduced distance for a given strength of the sphere-surface interaction. 
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The two branches define the distortion of the resonance peak as a function of d. U- gives 
the evolution of the resonance peak for frequency values below the resonance one and m+ 
for frequency values above the resonance. 

For the description of a Tapping experiment, the variation of the amplitude a as a function 
of the distance d is readily obtained by rewriting equ^l as follow : 



d 



± 



\ 



n 2/3 



3|g(l-M2)^^l/a2-«2| 



(17) 



Consequently, depending on the drive frequency and the drive amplitude (through the A, 
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dependance into bifurcations from a stable to a bistable state may occur leading to 
amplitude and phase jumps. From an experimental point of view, the conditions required for 
the appearance of the bifurcations have been extensively discussed in refs.— "2^. In particular, 
the use of drive frequencies lower than the resonance frequency favor the measurement of 
the bifurcations corresponding to the non-contact situations. 

In figsO and El are given the distortion of the resonance peak and the evolution of the 
amplitude as a function of d for an identical set of parameters. For large values of d, 
e.g. when the surface is far from the OTCS (point K), the non-linear effects are negligible 
and the peak keeps a well-defined Lorentzian shape (see equJKil with Ka = 0). When the 
OTCS is approached towards the surface, because the interaction is attractive, the resonance 
peak starts to distort towards the low frequencies. The distortion of the peak increases as 
d decreases. In the vicinity of the resonance, m < 1, a > 1 and for small values of d the 
branches u+ and m_ become very close. To mimic a Tapping experiment, the drive frequency 
is fixed to Udrive in figO With the parameters used, at di = 2 the oscillations properties are 
nearly identical to those at d infinite. When the peak starts to distort the amplitude and 
the phase (not shown) vary (L). In the present example, b chosen below the 

resonance, the amplitude first increases. When the OTCS is further approached it reaches 
an unstable branch (M) and jumps to the stable branch (N) such that the bistable structure 
of the oscillator can be experimentally observed. At a closer distance, the peak further 
distorts and the amplitude is reduced since it follows the variations of (point O on figO])- 
Then, when the oscillator is retracted, it follows the upper stable branch until the resonance 
value corresponds to the fixed frequency Udrive (point P on figEl) and then jumps down to 
the lower stable branch. The curve exhibits a hysteresis cycle (points M, N, P and L). On 



the figs|21 and El the branches that are supposed unstable are shown with dashed hues. 
D. Resonance frequency shift 

Using equHni the resonance frequency shift as a function of the distance d is obtained by 
setting a = 1. This former condition ensures the required condition for the NC-AFM mode. 
Thus, the normahzed frequency shift, {v — Vq) /z/q, is given by m — 1— : 



In figEl is given the frequency shift as a function of the distance d for the same set of 
parameters than the one of figsOandEl Following the previous discussion about the stability 
of the different parts of the resonance peak during the distortion, since the measure is 
performed as a function of d with a = 1, no bistable behavior can be observed. The 
OTCS follows the same branch m_ or m+ whose location is always stable whatever the peak 
distortion. 

Nevertheless, note that it should exist only one branch of variation for the frequency shift 
which is defined from the condition a = 1. But two branches are obtained as a consequence 
of the two solutions d^. When the peak is distorted, the branches u_ and u+ become very 
close as Q becomes large (see for instance, figOwith d^ = 1.012). Therefore even with an 
oscillation amplitude kept constant, question rises about the ability of the OTCS to remain 
on the same branch. Qualitatively, one may expect that around a = 1, the branch m_ is 
unstable and is stable (see figl21)- If this is true, any small fluctuation of the oscillation 
amplitude might produce a jump from one branch to the other one as discussed in ref.— . 
Since the branch m_ seems to be unstable, a jump to this branch should lead to an abrupt 
decrease of the amplitude, which in turn might produce an apparent abrupt decrease of the 
quality factor. Because such a jump should show accidents, both on the resonance frequency 
shift curve and on the damping signal, accidents which are, in most cases, not observed, it 
becomes useful to determine more accurately the stability of the two branches. 

The main aim of the variational method is to define a theoretical frame describing at the 
same time the evolution of all the variables : a, ip, u and d of the OTCS dynamics. Thus, 
it becomes useless to discuss about the stabihty of each kind of branch for the two dynamic 




(18) 
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modes since the stability of a given couple (a, or (a, u) as a function of d can be deduced 
for any other couple. The discussion will be made on the branches of the resonance peak. 

III. STABILITY CRITERIONS 

The stability of the branches u± of the resonance peak (see equ]16|) is obtained from 
equations of motion of the OTCS (see equJ14j). These equations are linearized around the 
stationary solution which will be now identified by the index "s" . At this stage, we fall into 
the well-known linear theory (see for instance refs.— and2^). Writing into a four dimensions 
matrix the linearized system corresponding to the generalized coordinates (a, ip) and their 
associated generalized velocities (a, 0), "we extract the eigenvalues and discuss the stability 
as a function of the sign, negative (stable solution) or positive (unstable solution), of their 
real part. The calculations are detailed in the appendix. 

The stability criterion for each one of the two branches is given by the following inequation 
(see equEHl) : 

A is proportional to (see appendix, equl22|). thus characterizes the non-linear attractive 
coupling term. Solving the associated equality gives the critical branch m^T** whose loca- 
tion regardless u± defines the stability of each branch. Unfortunately, m^!"** has no simple 
expression. Therefore the stability condition can not be exploited as is. Nevertheless a 
numerical routine allows to get its main features. The figElfay) shows the distortion of the 
resonance peak and the critical branch u^** numerically computed for each branch, u^** 
never crosses u+ and is always located below it. On the contrary, m^'* crosses U- twice and 
their relative position depends on the value of the amplitude which, in turn, is going to 
define two domains of stability. The figs\B(b) a.nd\^c) are zooms on the regions a and (3 of 
The intersection spots are exactly located where the curvature of m_ changes. Therefore 
it's worth discussing the stability as a function of the local curvature of the branches u±. 
and so introducing their derivative das/du±. As shown in the appendix, the inequation IT^ 
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can be summarized as follow 




and 



cos {(ps) > aj (2(5) 



< 



or 



(20) 



V 




and 



cos {(ps) < aj (2(5) 



II 



The aim of this former expression is to exhibit an explicit and particularly simple dependence 
of the stability of the branches as a function of their derivative : 

S-i' For the u+ branch, das/du+ being always negative and the associated value of the 
phase being always defined beyond —90° (see section ITm|l . thus cos (^9^) < 0, the criterion 
{ii) implies that the u+ branch is always stable, whatever the value of a^. 

S-^ Concerning m_, the sign of the derivative changes twice. For this branch, the phase is 
always defined above —90° which in turn means cos {(ps) > 0. Therefore on the lower part of 
the branch (small a), dag/du^ > and the criterion (i) indicates that the branch is locally 
stable. When dag/du- becomes negative (see Hg^b)), because cos{ips) is still positive, the 
criterion (i) is no more filled. As a consequence, U- is locally unstable and the instability 
is precisely located where the infinite tangent appears. On the upper part of the resonance 
peak, the curvature of u_ changes again and das/du_ > (see fig^c)), implying that it 
is again a locally stable domain. Thus the branch m_ exhibits two stable domains and one 
unstable. 

Note also that the resonance condition is deduced from dag/du = which implies 
cos{(fs) = CLs/ (2(5)- This equality is the usual resonance condition of a free harmonic 
oscillator. If = 1, e.g. without any coupling, the resonance phase is therefore 
(fis = arccos [1/ {2Q)]. For the OTCS we used, Q ^ 500, and so y?, ^ -90°. 

IV. DISCUSSION 

The previous criterions allow to conclude to the stability of the OTCS for each dynamic 
mode. It was shown that, for an attractive coupling, the branch was always stable and 
that the instability was controlled by u„. For frequencies lower than the resonance (branch 
M_), when the tangent dag/du- is positive, the branch is stable. Thus, there is a small 
domain close to the resonance value for which the U- branch remains stable. 
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For the Tapping, this result imphes that the OTCS is locked on a stable branch until 
das/du- ±00 (see figsO and 01 point M) which makes the amplitude jumping up to the 
upper stable branch u+ (point N on figOI) during the approach and jumping down to the 
lower stable one u_ (point P) during the retract. 

For the NC-AFM, the result of the present work shows that is always stable but that 
also a small domain of U- around the resonance value remains stable. If the resonance 
value would have been located at the point where das/du^ is infinite, an infinitely small 
fluctuation would have been able to generate a catastrophic behavior like large variations of 
the oscillation amplitude and lead to an abrupt increase of the damping signal as previously 
discussed and suggested in ref.— . 

Nevertheless, the size of the m_ stable domain is Q dependant. The more the Q factor 
is large, the more the size of the domain is reduced. The figEl illustrates the reduction of 
the size of the domain for Q = 500 and Q = 5000. With Q = 5000, the size of the domain 
is so weak that it nearly no appears on the scale of the figure. As a consequence, even if 
the previous discussion may help to understand why abrupt increases of the damping signal 
do not systematically occur, question remains unclear for the very high Q factors that can 
be obtained in ultra-high vacuum {Q > 10000). At this step, it's worth giving orders of 
magnitude. During a NC-AFM experiment, an electronic feedback loop keeps constant the 
amplitude of the OTCS and locks its phase at —90°. Therefore question rises about the 
size of the stable domain in phase around —90°. If any fluctuation around the locked value 
goes beyond the stable domain, the OTCS behavior becomes unstable. For Q = 500, the 
size of the stable domain is of about 1.5° whereas it's reduced to 0.15° for Q = 5000 (data 
not shown). Thus, if the electronic loop is able to control the phase locking with a better 
accuracy than 0.15°, the OTCS will be locked within a stable domain and in turn won't 
give rise to instabilities. In addition, since the ability of the electronic loop to control the 
oscillating behavior depends on the value of the quality factor—, the reduction of the domain 
of stability might not be a key parameter. 

Practically, during our experiments, even with quality factors larger than 10000, drastic 
variations of the oscillation amplitude are never observed. Thus, the main aim of the present 
work is to show that, if the oscillator is properly locked at the —90° value throughout an 
experiment, this value corresponds to a stable domain. 
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V. CONCLUSION 



This paper was a theoretical investigation of the stabihty of the non-hnear behavior of 
an oscillating tip-cantilever system close to a surface. A variational principle has allowed 
to get the temporal dependance of equations of motion of the oscillator as a function of 
the non-linear attractive coupling. The interaction potential chosen is a disperse Van der 
Waals one, calculated between a sphere and a plane. The stationary state is obtained and 
can be interpreted either in the Tapping mode or in the NC-AFM mode. The stability of 
the stationary state is analyzed in terms of distortion of the resonance peak as a function of 
the coupling. It is found that stability criterions can be expressed from a simple inequality 
involving the sign of the derivative of the curve. The branch associated to the frequencies 
larger than the resonance is always stable whereas the branch associated to the frequencies 
smaller than the resonance exhibits two stable domains. The instability appears when the 
branch exhibits an infinite tangent. This feature allows to re-interpret the instabilities 
appearing in Tapping mode and may help to understand why the NC-AFM mode is stable 
most of time. 

Appendix: Computation of the stability of the branches 

Let's note a = as {1 + ^) and ip = ips + p with ^ 1. The index s is attributed to the 
stationary solution. Keeping the terms of first order in equ lT^ equations of motion of the 
variations C, and p regardless and respectively, may be written as : 




(21) 



with : 




(22) 
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The system is solved by setting S = ^, and \1/ = p leading to a linear system of the fourth 
order : 



with : 



P 



P 



(23) 



/ 





1 \ 

1 

u/Q -1/Q 2u 



(24) 



\ -u/Q - cos (ifis) / (QcLs) -2m -l/Q J 
The eigenvalues of the matrix are obtained by solving the characteristic polynom given by 
P = det (M — XT). P can then be written as : 



P={X^ + X/Q + M) (A" + X/Q + N) 



By identification : 



M + AT = + 1 - A 



cos (v^s) 



MN = {u/Qf - (m2 - 1 + A) 



cos [ips] 

Qa-s 



The characteristic equation P = is then equivalent to the following system 



+ A/Q + M = 
X^ + X/Q + N = Q 



Ai,: 



^3,4 



-1/g ± sJ{l/Qf-AM^ 12 
-l/Q±,J{l/Qf-m^ /2 



The stable solutions are the ones given by 3? (Aj) < 0^, thus 



(25) 



(26) 



(27) 



IIQ> \I{IIQY -AM 

and M > and iV > 



l/Q>\Kl/QY-AN 



(28) 
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According to relationships 123 two conditions are necessary to fill equ|2Hl and in turn to 
ensure the stability of the solutions : MN > and M + N > 0. 

S-i>Let's first consider MN > which is the main of the both (see below) : 



MiV=f^V-(n^-l + A)^^>0 (29) 

The equation MN = can be numerically solved using a Maple routine. Nevertheless a 
tractable stability criterion requires to write in a different way the expression [201 Using the 



relationship cos (yj^) = iyl — (uas)'^ (see equlT^. the two coupled equations of the sine 
and cosine of the phase of the stationary state imply : 



with : 



Therefore 



G {as,u) = Qtts {l - u^) - g{as) T - (uas) =0, (30) 



dG (os, u) = da^G {as, u) das + duG {as, u) du = 0, (32) 



and so : 



das duG {as, u) 



du da,G{as,u) 

The calculations lead to : 



QJ Qa_ 
The denominator is exactly the product MN , thus : 



du 



(33) 



das_ _ u_ 2cos{ips) - as/Q . . 



MiV>0 4»^x '°°''''^:'-''""^ >0 (35) 

Q das 
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u/Q being always positive, the stability condition MN > is reduced to : 



2cos(v9s) - ds/Q 




and 



cos {ips) > as/ (2(3) 



dag 
du 



> ^ < 



or 



(36) 




and 



cos(v?s) < as/ (2(5) 



S->Let's now consider the condition M + N > : 



M + N = 3u^ + l-A + 



cos (ips) 

Qas 



> 



(37) 



Considering that varies within the range [0..1 + e], with 1 ^ e > 0, it's straightforward 
to show that the inequality is always filled so that the stability criterions are only given by 
the condition MN > and equHHl 
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FIG. 1: Sketch of the OTCS and notations used. 
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FIG. 2: Evolution of the resonance peak computed from equation (16) for three values of the 
distance, di = 2, d2 = 1-11 and = 1.012. The numerical parameters are = 20 nm, Q = 400 
and Ka = 8 X 10^'*. For an attractive coupling, the peak is more and more distorted towards the 
low frequencies as d is reduced, e.g. the surface is approached and lead to the bifurcations observed 
on the tapping curve (Fig. 3, point M) when the drive frequency, Udrive = 0.9985, is chosen below 
the resonance and to the variations of the resonance frequency shift (Fig. 4). For each value of d, 
the unstable domains of n_ are shown with dashed lines. 
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FIG. 3: Variations of the amplitude as a function of the distance, e.g. approach-retract curve in 
the tapping mode computed from equation (17). The numerical parameters are the same than in 
Figure 1 and the drive frequency chosen is Udrive = 0.9985. The curve exhibits a hysteresis cycle 
(MNPL) due to the non-linear coupling that characterizes bifurcations (points M and P) from a 
monostable to a bistable state (see text). The stable domains of the branches are shown with 
continuous lines and the unstable domain with dashed line. 
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FIG. 4: Variations of the frequency shift of the resonance peak of the oscillator as a function of 
the distance, e.g. approach-retract curve in NC-AFM mode computed from equation (18). The 
numerical parameters are the same than in Figure 1. It's predicted that the curve is stable with d. 
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FIG. 5: (a) Distortion of the resonance peak computed from equation (16). The numerical param- 
eters are d = 1.05, Aq = 10 nm, Q = 500 and = 2.5 x 10^'*. The critical branches u^T** are 
calculated from equation (19). n+ is always located above its critical branch u^** and in turn is 
always stable, u^** crosses U- twice. This leads to define three domains to describe the stability 
of the branch. The domains are defined between the spots where the derivative da/du-^ diverges 
(see text), (b) Zoom in the region a of n_. The stability criterion foresee that below u"^**, n_ is 
stable and unstable above. This is illustrated by the dashed lines, (c) Zoom in the region (3 of n_. 
As da/du- diverges again, it defines a new domain of U- which is predicted to be stable. 
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FIG. 6: Zoom on the same scale than Figure 5c of the distortion of the resonance peak for two 
values of the quality factor of the OTCS, Q = 500 and Q = 5000. The numerical parameters are 
the same than in Figures 5. The unstable domain of u_ is still shown with dashed lines but the 
size of the second stable domain is drastically reduced for the larger value of Q so that it nearly no 
appears on this scale. The associated phase variations regardless —90° (not shown) are of about 
1.5° and 0.15° for Q = 500 and 5000, respectively (see text). 
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